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Increasing the spin imbalance in superconductors can spatially modulate the gap by forming 
Cooper pairs with finite momentum. For large imbalances compared to the Fermi energy, the 
inhomogeneous FFLO superconductor ultimately becomes a normal metal. There is mounting 
experimental evidence for this scenario in 2D organic superconductors in large in-plane magnetic 
fields; this is complemented by ongoing efforts to realize this scenario in coupled tubes of atomic 
Fermi gases with spin imbalance. Yet, a theory for the phase transition from a metal to an FFLO 
superconductor has not been developed so far and the universality class has remained unknown. 
Here we propose and analyze a spin imbalance driven quantum critical point between a 2D metal 
and an FFLO phase in anisotropic electron systems. We derive the effective action for electrons 
and bosonic FFLO pairs at this quantum phase transition. Using this action, we predict non-Fermi 
liquid behavior and the absence of quasi-particles at a discrete set of hot spots on the Fermi surfaces. 
This results in strange power-laws in thermodynamics and response functions, which are testable 
with existing experimental set-ups on 2D organic superconductors and may also serve as signatures 
of the elusive FFLO phase itself. The proposed universality class is distinct from previously known 
quantum critical metals and, because its critical fluctuations appear already in the pairing channel, 
a promising candidate for naked metallic quantum criticality over extended temperature ranges. 

PACS numbers: 74.40.Kb, 71.10.Hf, 05.30.Rt, 74.70.Kn 


I. INTRODUCTION 


The concept of electronic “quasiparticles” moving 
through a crystal and scattering off each other at rates 
smaller than their typical kinetic energy has -together 
with fermionic quantum statistics- led to a satisfactory 
understanding of the electrical conductance and the for¬ 
mation of superconductivity observed in many metals 
at low temperatures. However, an increasing number 
of observations in technologically relevant compounds 
indicate intriguing “strange metal” phases^, with long- 
ranged quantum entanglement^, where interactions de¬ 
stroy the electronic quasiparticles, and where electrical 
currents flow at anomalously slow rates^. 

Here, we propose a new strange metal phase associated 
with an underlying quantum critical point (red region in 
Fig. 1 ) in anisotropic electronic systems at the onset of 
inhomogeneous FFLO (Fulde-Ferrell-Larkin-Ovchinikov) 
superconductivity^^® in two spatial dimensions. We ar¬ 
gue that organic superconductors'®’^^ are promising can¬ 
didates hosting this new phase of matter in which elec¬ 
tronic quasiparticles are destroyed over regions of the 
Fermi surface from scattering off FFLO waves carrying 
a hnite momentum. As a consequence, specific heat and 
NMR (nuclear magnetic resonance) relaxation rates be¬ 
come non-Fermi liquid power-laws, whose exponents we 
compute. We propose that a new round of data-taking 
with existing experimental set-ups'^’® on for example the 
compound k—(BEDT-TTF) 2Cu(NCS)2 could provide a 
near-future experimental test of our predictions. 



FIG. 1: New strange metal phase with non-Fermi liquid be¬ 
havior (red region) at the quantum critical point for the on¬ 
set of inhomogeneous FFLO superconductivity in anisotropic 
electronic systems in an in-plane magnetic field (h) in two 
spatial dimensions. In the organic conductor k—(BEDT- 
TTF) 2 Cu(NCS) 2 , recent evidence places the onset of FFLO 
correlations (green region) to high in-plane magnetic fields 
h ~ 24 — SOT and low temperatures (T < 4K)'*“®. The black- 
dashed line is crossover scale Tcross explained in the text. 

A. FFLO phase in organic superconductors 

Organic superconductors have emerged as leading can¬ 
didates for observing FFLO superconducting states in 
electronic systems^“®’^®“^®; complementary efforts are un¬ 
dertaken in heavy Fermion^® and iron-based^^’^® com¬ 
pounds, as well as in ultracold atoms^®^®^. The crystals 
of organic superconductors are grown relatively clean. 
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with a mean-free path larger than the coherence length 
of the pairs carrying a hnite momentum Qfflo- For 
in-plane magnetic fields, superconductivity can be de¬ 
stroyed by a Zeeman splitting before orbital effects be¬ 
come relevant. For an isotropic superconductor with no 
orbital effects, the critical Pauli limit field Hp is reached 
when the Zeeman energy ^pHp > Aly/^ overcomes the 
gap^^’^^. A Zeeman driven transition^"^ requires an or¬ 
bital critical field F/^orb larger than Hp. The Maki pa¬ 
rameters a = V^Fforb/F^p must thus be larger than one, 
as in /c-(BEDT-TTF) 2 Cu(NCS )2 where a ~ SA 

Due to only very weak and/or incoherent elec¬ 
tron motion between the two-dimensional layers 
(^p/^interlayer ^ 3700 for k-(BEDT-TTF) 2 Cu(NCS) 2 '‘ 
and ^ 500 in the Bechgaard salt 

(TMTSF) 2 C 102 ^^’^^), out-of-plane orbital currents dis¬ 
sipate quickly. This is essentially true for the two 
compounds mentioned previously and also in the 
/3"-phase compound / 3 "-(BEDT-TTF) 2 SF 5 CH 2 CF 2 S 03 , 
whose phase diagram and NMR spectra do not depend 
on the in-plane orientation of the B-field®^. 

Fig. 2 shows three approximate Fermi surfaces from 
FFLO-candidate organic superconductors that share flat 
regions to which Qfflo would couple to preferentially. 
Qfflo are also incommensurate with the underlying 
crystal momenta. We expect our main results to be (i) 
universal across the various Fermi surface topologies, as 
well as (ii) relevant for both uni-directional (discussed 
later) and bi-directional FFLO modulation provided that 
the Fermi surface has (at least locally) nested parts which 
can -and will- be connected by fluctuating Qfflo order. 

II. MODEL 

To understand the main features of prototypical 
anisotropic 2D Fermi surfaces depicted in Fig. 2, let us 
consider a simplified model with a spin-dependent single 
particle dispersion 

^o-(k) = — 2ty cos{dky) — fJ. — crh. (1) 

Here, m is the effective mass for the motion along the 
chains in x-direction. The chains are coupled by a weak 
hopping matrix element ty as is the case for example in 
the Bechgaard salt (TMTSF)2C102^^’^^. There, the hop¬ 
ping parallel to chains is tx ^ 1340 K, between the chains 
ty ^ 134 K and much weaker in the out-of-plane direction 
tz ^ 2.6 the latter is therefore neglected in ^cr(k). 
Orbital effects from t;,-hopping influence the (material- 
specific and non-universal) location of For the 

K-compound the interlayer motion is very small and for 
the /^''-compounds even incoherent, and plays therefore 
no role^®. d is the distance between the chains (^ 7A in 
(TMTSF) 2 C 102 )), ^ is the chemical potential, and h the 
Zeeman field which splits the energies for spin-up {a = 1) 
and spin-down electrons (cr = —1). The key feature of 

(k) is a nested pair of so called hot spots (indicated by 
black points in Fig. 2) connected by Qfflo; tke universal 


“Qfflo Qfflo 



FIG. 2: Approximate Fermi surfaces geometry of candidate 
materials for the FFLO phase and associated strange metal, 
(a): K— (BEDT-TTF) 2 Cu(NCS )2 with closed elliptic hole 
pockets corresponding to a orbit frequency of about 600 Tesla'^ 
and open parts roughly aligned with the c-axis^®. Therefore 
an in-plane magnetic field of F = 30 T results in a ~ 5% 
mismatch (exaggerated in the figure) of the spin-up and spin- 
down Fermi surfaces. Preferential FFLO ordering possibilities 
include a uni-directional modulation in c-direction along the 
short axis of the hole pocket, where Qfflo would nest and 
connect the (Zeeman-split) sides of the hole pocket, (b): /?”" 
(BEDT-TTF) 2 SF 5 CH 2 CF 2 S 03 with tilted hole pockets. A 
preferential two-dimensional Qfflo, with components along 
the b and c axis, would again connect the weakly curved parts 
of the hole pocket at a right angle. In both the k- and /?”- 
phase interlayer transport is weak and incoherent such that 
orbital effects can be neglected, (c): Open Zeeman-split 
Fermi sheets characteristic for the two-dimensional arrays of 
conducting chains of the Bechgaard salt (TMTSF)2C102^®’^"^ 
in an in-plane magnetic field. The precise warping of the 
Fermi sheets also depends on the applied pressure^^. 


feature it shares with the Fermi surface topologies also of 
k-(BEDT-TTF) 2 Cu(NCS )2 for short-axis FFLO mod¬ 
ulation, and / 3 "-(BEDT-TTF) 2 SF 5 CH 2 CF 2 S 03 with di¬ 
agonal two-dimensional FFLO modulation as per Fig. 2. 
The most relevant interaction for the onset of spin-singlet 
FFLO pairing is a short-range attraction g: 

Hint =-9 J 

where the electrons are represented by anti-commuting 
operators "(/J., ipa- For our calculations below, we can 
assume a weak coupling scenario, where the microscopic 
origin of the pairing interaction plays no role®^®'^. The 
Fermi surface mismatch between spin-up f and spin-down 
/ electrons frustrates conventional homogeneous BCS 
pairing, where the Cooper pairs Abcs = (^-k.tV'kq) 
carry zero net momentum and the superconducting gap 
function is homogeneous in space. As analyzed in Fig. 3 
for the specific (k) given above, the mismatched Fermi 
surfaces force the electrons to pair at a finite momentum 
Qfflo = (QojO) leading, here in x-direction along the 
chains, to a modulated order Afflo = Aq cos(QFFLO’r), 
like at a liquid crystal transition^®’^®. 
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III. NOVEL EFFECTIVE ACTION 


At zero temperature T = 0, the FFLO state forms out 
of the imbalanced metal at large magnetic fields upon de¬ 
creasing h via a continuous quantum phase transition at 
h = /iQCP in Fig. I (see also Appendix). Gaps open at the 
hot spots, i.e. at the discrete points on the Fermi surface 
connected by Qfflo- Our aim is now to describe quan¬ 
tum fluctuations of electrons around the hot spots that 
can scatter off the incipient FFLO order with momentum 
transfer ±QFFLO + <5q between hot spots. We parameter¬ 
ize the electron operator by fermionic quantum fields 
k) corresponding to the right and left part of the 
Fermi surface, as indicated in Fig. 2 (c). The construc¬ 
tion goes completely analogously for the Fermi surfaces 
Fig. 2 (a,b). Fluctuations of the FFLO order parame¬ 
ter are described by one bosonic complex quantum field 
A(r, k)i 2 per hot spot pair, with r the (imaginary) time 
coordinate. Fluctuations of both, the electrons and the 
FFLO field, are encoded in the partition function 

Z = J exp(-5) , (3) 

where the exponential of the action S = dr f d^r C 
weights the field configurations according to: 


^ = 5 1 ^* 1 ^ + 51 ivid^ + 


2 = 1,2 


j^R,L 


2m., 




- 9 V'f) + h.c] 


( 4 ) 


Here we have expanded the electronic dispersion around 
the hot spots ^^’^(k) = + ky/2my, with x,y be¬ 

ing the direction orthogonal and parallel to the Fermi 
surface, respectively, my is (inversely) proportional to 
the curvature of the Fermi surface (and thus to ty) and 
related to the applied pressure in the Bechgaard salt 
(TMTSF)2C102^^’^'^’^^. Further, we note that in the ab¬ 
sence of time-reversal symmetry Eq. (4) is not invariant 
under the spin-flip transformation The fermion- 

boson vertex g converts bosons with momentum Qfflo 
into electron pairs, one at each hot spot of a given pair. 
The finite curvatures of the Fermi surface allow to re¬ 
strict Eq. (4) to the hot spots provided temperatures are 
smaller than Tcmss ^ tyh/ep (cf. Fig. I). 


A. Distinct universality class 

The electronic propagators in Eq. (4) bear some sim¬ 
ilarity to those in two-patch theories for metals cou¬ 
pled to nematic and C/(l)-gauge field fluctuations^°’^^ 
as well as for the onset of incommensurate 2k-p 
antiferromagnetism^^. An important difference to these 
latter cases is that the electronic interactions in the 
particle-hole channel and that they preserve time-reversal 



FIG. 3: (a): Pairing susceptibility n(tj„ = 0, Q) at zero 
frequency and finite momentum for the Fermi surfaces c) in 
Fig. 2. The position of the maximal peaks at Qfflo = (Qo, 0) 
determine the, here unidirectional, modulation wave vector 
of the FFLO-phase. (b): cuts of the pairing susceptibility 
around the peak at Qfflo, displaying the asymmetrically- 
sloped ridges along characteristic of incommensurate 
symmetry-breaking. To capture quantum fluctuations around 
the peaks ±Qfflo + <5q, we later parameterize the pair¬ 
ing fluctuation propagator with a linear slope from the 
“outside” and a non-anal 3 rtic square-root from the “inside” 
ridge. Along Qy the peak is symmetric. The susceptibil- 

ity n(u;„,q) = ^ Ek - e^(k)+e^(q-k)-i.,„ - involves a 

sum over available states in momentum space (numer¬ 
ator) weighted by the inverse of their energies (denominator). 
Here npiE) = (1 -|- exp(/?E))“^ is the Fermi function and 
P = llksT the inverse temperature. Parameter values for the 
plot are ty = 0.5 fmd^, g = 3.3 tP/mdP, h = 1.0 fmd^. 


symmetry. In contrast, Eq. (4) defines a new universal¬ 
ity class as it describes singular electronic interactions in 
the (FFLO) pairing channel and the Zeeman field breaks 
time-reversal symmetry. 


IV. ANOMALOUS QUANTUM DYNAMICS 


The first step toward understanding quantum fluctua¬ 
tions is to compute the dynamics of the bosonic FFLO 
pairing field I?i=i(T, r) = (Ai(t, r)A|(0, 0)), which it 
picks up from emission and absorption of electron pairs. 
In frequency and momentum space, the interest lies in 
the slopes and direction dependence of the pairing sus¬ 
ceptibility close to the maxima in Fig. 3. The leading 
contribution to I?i(a;„,k) is given by the one-loop bub¬ 
ble diagram evaluated in the Appendix. For the criti¬ 
cal properties of interest here, it suffices to expand the 
(massless) fluctuation propagator directly at the quan¬ 
tum critical point for small frequencies a;„ and momenta 
TQfflo + ^k. The result in the physically relevant small 
imbalance limit, Sv = — vi ^ v = (i;^ -I- u^)/2 (re¬ 

call that an in-plane magnetic field of 30 T leads to a 
~ 5% mismatch in spin-up and down Fermi surfaces in 
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k-(BEDT-TTF) 2 Cu(NCS) 2 '‘) is: 


(w„,k) = 


Attv\5v/v\ 


ky Sv 

2Re-i / — ——h dvkx H- iu}„ 


2my 


kl 6v 

+B- -h C — 

ZniyV V 


( 5 ) 


I. We note this value is also obtained in nematic metals 
at the one-loop level'^° but receives quantitative correc¬ 
tions at higher loops^^’^^. We expect diagrammatic ex¬ 
pansions of the path integral of Eq. (4) to be different 
to the nematic metal, however. One-loop corrections to 
the fermion-boson vertex, for example, vanish here in the 
EELO case. 


The key feature here is the square root, introducing a 
non-analyticity; the usual poles of coherent excitations 
are substituted by a branch cut in the complex plane. 
Physically, this reflects the overdamped, dissipative dy¬ 
namics of the order parameter due to decay into electron 
pairs. The subleading corrections to the square root sin¬ 
gularity proportional to are the first analytic cor¬ 

rections which capture the asymmetry in the slopes of 
the peaks in Fig. 3. Eq. (5) bears some resemblance to 
the one-loop results for the order parameter dynamics in 
the incommensurate antiferromagnetic and charge order 
case^^’^^, which are singular fluctuations in the particle- 
hole channel of different physical origin. 

A. Non-Fermi liquid behavior of electrons 

The peculiar dynamics Eq. (5) of the EELO pairing 
fluctuations strongly back-acts on the electrons close to 
the hot spots. This renormalizes the electron propagator, 
say, at the L f hotspot, GL^{T,r) = {tpL^{T,r)tpl^^{0,0)) 

whose Fourier transform reads k)] ^ — 

^;^(k) The inverse quasi-particle lifetime 

T”^ at the hot spots for propagators of the form (5) can 
be computed analytically^^. At the QCP, evaluation of 
the one-loop self-energy with the dynamical pairing prop¬ 
agator yields (see Appendix) 

exhibits non-Fermi liquid behavior with a fractional ex¬ 
ponent 2/3 due to the interactions with the singular 
EELO pairing fluctuations (recall that in a Fermi liq- 
uid Pp/ ~ w^). The strength of the anomalous decay in 
Eq. (6) is proportional to the difference Sv of the Fermi 
velocities which itself is proportional to the applied Zee- 
man field. In the language of critical phenomena Eq. (6) 
implies an anomalous frequency dimension tJt = 1/3 as 
the bare propagator is linear in frequency^^^^®. Probably 
weaker non-Fermi liquid behavior may also appear in ho¬ 
mogeneous Fermi systems with Zeeman-imbalance^^^^®. 
A similar computation for the electron self-energy in kx 
and ky direction yields a much weaker renormalization 
such that the corresponding ry^’s vanish at the one-loop 
level. For direct loop computations (in absence of in¬ 
frared cutoffs), the fermionic dynamic exponent is defined 
as = \Zy*‘ Slid deviates from the Fermi liquid value 

(zpL = 1) and instead attains the value z^p^Q one-ioop ~ 


V. PREDICTIONS FOR NEW EXPERIMENTS 

The general scaling form of the electronic specific heat 
in the quantum critical “fan” in Fig. 1 is^^’^^ 

CejT ~ , (7) 

where the last equality has used ^pppo one-ioop = 
d = 2, and 0 = 1 in case hyperscaling is violated. Recall 
that in a Fermi liquid CejT is independent of tempera¬ 
ture. For K—(BEDT-TTF) 2 Cu(NCS) 2 , the specific heat 
has been measured by Lortz et over a range of B- 

fields including values towards the upper phase-boundary 
of the FFLO phase: an upswing with decreasing T incom¬ 
patible with Fermi liquid behavior can be seen at ft. = 28 
T but new rounds of data taking at high magnetic fields 
are needed to resolve this. 

A more recent experiment, where signatures of FFLO 
order have been detected in the compound k—(BFDT- 
TTF) 2 Cu(NCS )2 was performed by Mayaffre et al. who 
have measured the NMR relaxation rate I/(rir) at high 
magnetic fields including data set at ft = 22, 27 Tesla 
(see Fig. 2 of Ref. 6). They associate the observed up¬ 
swing at ft = 22 Tesla with bound states of quasiparticles 
at the nodes of an FFLO order parameter. At ft = 27 
Tesla, however, I/(TiT) seems to decrease for tempera¬ 
tures smaller than 3 Kelvin. Our strange metal sccenario 
predicts a power law decrease of the form I/(TiT) ~ 
which is suppressed compared to a constant Fermi liquid 
relaxation rate. New rounds of data taking over extended 
temperature and higher magnetic fields should resolve 
this. The power-law exponent in I/(TiT) is determined 
by the scaling behavior of the density of states 

^ ~ 11 du;N^°\u;)N^°\u;)nAu^) [I - ( 8 ) 

with = f For the critical elec¬ 

trons at the hot spots the associated spectral function is 
expected to fullfill lu/T scaling'^^ds 



where the anomalous scaling in spatial directions x and 
y can take different finite values but do not influence Zf 
at the one-loop level. The same temporal anomalous di¬ 
mension 77 t- = 1/3 that appears in the self-energy Eq. (6) 
leads to 7V^°*(w) ~ (see Appendix) and thus to the 
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exponent 2/3 in Eq. ( 8 ) for electrons from the immediate 
vicinity of the hot spots. As in the nematic metal, we ex¬ 
pect higher-loop corrections to further suppress the den¬ 
sity of states^^. In contrast to the nematic metal, how¬ 
ever, where all points on the Fermi surface can be taken 
to be equivalent^^’^®, here in the FFLO strange metal a 
discrete pair of hot spots is singled out. To account for 
electrons from “colder” parts of the Fermi surface one 
should add a constant Fermi liquid piece to the density 
of states Nc,{uj) = -I- iV^°‘(w) and NMR relaxation 

rate when comparing with experiments. A proper anal¬ 
ysis of crossovers and potential contributions violating 
uj/T scaling^®d7 j.equjj .0 momentum-resolved renormal¬ 
ization group techniques. 

VI. SUMMARY 

The FFLO strange metal phase proposed here pro¬ 
vides a new universality class for the study of critical 
thermodynamics and transport without electronic quasi¬ 


particles. Our results open up the possibility to targeted 
new rounds of data taking in the organic superconduc¬ 
tors and may facilitate access to a naked quantum criti¬ 
cal metal with critical fluctuations in the pairing channel 
over extended temperature ranges. 
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Microscopically, the attractive interactions may thus be 
either due to a phonon-exchange mechanism or due to spin- 
fluctuation mediated pairing. 

Note also that a model with the simple single particle dis¬ 
persion Co-(k) of Fig. 2(c) and an adjustable attraction can 
be realized with ultracold atoms, where both the band 
structure and strength of the attraction can be engineered 
via optical lattices and Feshbach resonances^®’^^’'^®^®^. 
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Appendix A: Continuous mean-field FFLO transition 


We decouple the interaction term (Eq. (2) of the main text) in the pairing channel and write the partition function 
as the following functional integral 


Z = 




(Al) 


for the complex Grassmann field if) and the bosonic field A, with the euclidean interaction action 


^int = 


dr I d^r j^i|A(r,r)p - (A*(t, r)'!/';(r,r)'!/'t('r,r) + A(r,r)'(/>^(r,r)V’;(r,r)) 


An effective Ginzburg-Landau theory for the pairing field is obtained by integrating out the fermionic fields. Written 
in Fourier space the latter reads 


|A(w„,k)|2 - Trln/3G-i , 

9 . , 1 , 


(A2) 


where Tr = tr, with tr meaning the trace in the 2x2 spinor space, the bosonic(fermionic) Mastubara frequency 

LOn = 27rn//3(^'„ = 7 r( 2 n -|- l)//3), /3 the inverse temperature, and with the matrix propagator 


G-G 


,k') 


- CT(k)) A(:^„ - - k') \ 

- J^n,k'- k) +C;(-k)) y ’ 


(A3) 


where ^^(k) is the fermionic dispersion relation. Now, upon expanding the logarithm in powers of A, we can compute 
the coefficients of the effective Ginzburg-Landau theory. In particular, the second order term contains the pairing 
susceptibility 


n(a;„,q) = 

k 


1 - nF(/3G(q - k)) - nF(/3G(k)) 
G(k) -bG(q- k) - iuJn 


(A4) 


Specifying now to the mean-field FFLO transition with favoured momentum |Qfflo| we can express the paring field 
as 

A(cu7.j, q) y Ao(0, Qfflo) [^4'Qfflo^u,Qfflo(^^) A (27r) j . (A5) 

The d's are the amplitudes of the FFLO-order parameter, that are chosen to be equal and real. 

The quadratic Landau coefficient is thus given by 

02 = 1/ff — n(u;„ = 0, Qfflo) • 

The quartic Landau 04 coefficient results from the expansion of the logarithm to fourth order in A and contains 
all possible contractions of four fermionic propagators with the four external momenta fixed at ± Qfflo zero 
frequency. It must be numerically calculated. All contractions diverge when T —)■ 0, indicating the need to devise a 
renormalization group procedure for the problem at very low temperatures. For generic choices of parameters 

lim 04 > 0 (A 6 ) 

T->0 

indicative of a continuous phase transition at least at the mean-field level. 


Appendix B: Pairing fluctuation propagator, fermion propagator and density of states 


In this section we compute the one-loop corrections to the pairing and fermion propagators, leading to the Eqs. (5) 
and ( 6 ) of the main text. We use the low-energy Lagrangian from Eq. (4) of the main text, whose Fermi surface 
structure and one-loop diagrams are shown in Fig. 4. 

The inverse pairing fluctuation propagator diagrammatically shown in the upper row of Fig. 4(b) contains the bare 
part, simply given by the Landau coefficient 02 (vanishing at the critical point), plus the pairing susceptibility: 


Ili —1 (cu,^, q) — 


4 51 k)G^°^^(w„ - q - k) 


1 1 1 

^ ^ - Cf (k) i{uJn - Vn) “ ^^(q “ k) 


(Bl) 




FIG. 4: Elements of the low-energy quantum field theory for a single hot-spot of the Lagrangian given in Eq.(3) of the main 
text, (a): Fermi surfaces for the two species with dispersion relation + ky/2my. A finite curvature is present 

and the Fermi velocities of the two species are different, (b): Feynman diagrams at one-loop order. Upper row: propagator 
for the pairing fluctuations, containing the susceptibility given in Eq. (Bl). Lower row: self-energy correction to the electron 
dispersion and quasi-particle lifetime, given in Eq. (B3). 


We first perform the Matsubara sum in Eq. (Bl). By restricting to T = 0, we integrate first over kx and then over 
ky, always discarding the unphysical UV-divergent terms, and get 


Bj—i(ca^,k) — 




TO,, 


27r(r;^ -|- vi) 



- —+ (— - i){v^kx +iujn)+ 

vi 2my vi 



-p (^ - l)(^-y^kx + iUJn) 

2my 


(B2) 


which in the small imbalance limit <C u = (t;^ -|- vi)/2 gives Eq. (5) of the main text since k) = 

—11^=1 (w„,k) at the critical point. In the main text the first analytical corrections coming from higher order terms 
in the fermionic dispersion have been included in 11^=1 (w„,k). 

In order to compute the non-Fermi-liquid exponent we consider the quasi-particle decay rate given by the imaginary 
part of the analytically continued self-energy. At the one-loop level the fermionic self-energy is given by (Feynman 
diagram in the lower row of Fig. 4(b)) 


SLt(^'n,q) = 


^ ^ A=i(w„,k)G^](a;„ - z/„,k - q) . 


(B3) 


After analytical continuing —>■ w -|- zO^ and taking the imaginary part we get 

ImE)®j(a;, q) = - ^ [npi^Riik)) + nsi^Riik) + w)] ImH™‘(^i^;(k) -f w, k -h q) . (B4) 

k 


We now set T = 0, q = 0, and first integrate over kx and then over ky, to get Eq. (6) of the main text. By computing 
the self-energy (B3) at finite momentum we obtain no q-dependence at leading order, once we expand the pairing 
propagator (Eq. (5) in the main text) consistently with the small imbalance limit, to get the Landau-damped form 


\w„,k) 


V \/2my 
A'kv\5v/v\ 




-\5vlv\\ujn\ ^ ^ kl 


\ky 


2myV 


(B5) 


Within this approximate one-loop framework we have therefore no anomalous dimension for the fermion momenta: 
Vx = r]y = 0. 

By adding the self-energy correction from Eq. (6) in the main text to the retarded fermion propagator, we can 
compute the spectral function A^(w,k) = — ImG)(l|:(a;, k)/7r which satisfies the scaling relation of Eq. (9) in the main 
text with rjx = 1/3. The density of states 


N^{uj,T) 


(B6) 



























9 


is then obtained by integrating the spectral function over kx, ky after rescaling both coordinates in order to eliminate 
the w-dependence from the integrand. This procedure is based on the idea that we have to put a UV-cutoff on both 
the kx and the ky integrals since we are dealing with a low-energy theory in the proximity of the hot-spot. Accordingly, 
we add to the hot-spot contribution (B6) (called in the main text) a w-independent shift coming from the “cold” 
part of the Fermi surface (called in the main text), which indeed should behave like a Fermi liquid. However, 

in this way the exponent in Eq. (B6) can be strongly overestimated. 



